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INTRODUCTION 

Measurement  of  scattered  electromagnetic  field  and  further  permittivity  or  permeability 
reconstruction  based  on  experimental  data  and  adequate  mathematical  model  is  the  key 
issue  in  dielectric  materials  study  [1,2].  Here,  accuracy  of  measurements  and  adequacy 
of  mathematical  model  is  of  principal  importance.  Today,  special  attention  is  attracted 

to  the  study  of  thin  dielectric  films  with  10  3 +10  6 . In  this  paper  we  consider 

a resonator  that  can  be  used  for  a thin  film  study,  its  electromagnetic  model,  and 
advantages  and  capacities  of  corresponding  numerical  algorithm. 

DESCRIPTION  OF  THE  ALGORITHM 

Our  solution  to  inverse  problem  is  based  on  accurate 
and  efficient  solution  to  direct  problem.  Suppose  that 
the  diffraction  of  one  of  the  modes  of  small  circular 
waveguide  by  the  chamber  (Fig.l)  is  reduced  to  a 
functional  relation  Y = f{x).  Here  X = (x,.x:,x3....xm) 
is  the  set  of  input  data  (frequency,  S,  geometry, 
amplitude  of  incident  wave  A , etc.)  and 

Y - (y;,_y,,...,y„)  is  the  set  of  output  data  - reflection 

coefficients  RrR Ri; , normalized  by  A.  Assume  that 

a part  of  input  data  is  known  (frequency  and  geometry) 
and  given  by  the  values  x] , x2 , . . . x"  (/  < in)  from  the 

set  X of  all  input  variables.  Suppose  also  that  values 
Y"  - (y'l, y'i.... y°„)  of  output  variables  are  known.  The 
problem  of  the  model  identification  is  reduced  to  a 
necessity  of  finding  x,+j xm  from  equation 

r *„)  0) 

In  our  case,  we  have  to  find  unknown  dielectric  constants  = s'  + is" . The  accuracy  of 
the  parameter  evaluation  depends  on  several  factors: 

1.  The  error  of  reproduction  of  the  relation  between  output  data  Y,  known  input  data 
Xj,Xj,Xj...,xm,  and  unknown  parameters  in  the  form  of  equation  (1),  that  is  called 
inadequacy  of  the  model  to  the  phenomena  simulated. 

2.  Errors  of  measuring  the  known  parameter  values,  that  is  Y",  x" . a\  ....  x, . 

3.  Errors  of  the  numerical  algorithm  applied  to  solving  equation  (1 ). 
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As  a functional  relation  we  choose  the  solution  of  the  considered  problem  obtained  by 
the  semi-inversion  method  [4,  6].  Here,  we  reduce  it  to  numerical  solution  of  a matrix 
equation  of  the  second  kind  rm  - HRM  = B with  exponential  convergence.  This  approach 
allows  us  to  find  out  with  given  accuracy  the  amplitudes  of  reflected  modes 
Rf  ,r% if  the  frequency,  geometry,  and  dielectric  constant  a = s'  + is"  are  known, 

so  that  R"  = fJ/c,s,A)  depends  on  geometry,  frequency  parameter  k = ajX  (a  is 
characteristic  dimension,  see  Fig.  1),  and  e . The  accuracy  of  calculation  of  Rn  is  limited 

by  the  capacity  of  the  computer  utilized.  For  the  solution  of  direct  and  inverse  problems 
we  used  the  ideas  of  [7,  3-5].  We  pose  the  inverse  problem  as  a minimization  one,  see 
[6,3,4],  After  caring  out  numerical  investigation  of  various  types  of  functionals 
according  the  scheme  of  [4],  we  arrive  at  the  conclusion  that  the  most  efficient  is 

«W=iz  (2) 

m-l  n=I 

where  F(w)  = |wj  . Due  to  analyticity  of  functions  R^{k,s,A)  in  k and  £ , the 

functional  in  (2)  (within  the  considered  level  of  input  errors)  does  not  require 
regularization  and  we  can  put  (2)  a = 0 . 

NUMERICAL  EXPERIMENTS 

In  the  parameter  reconstruction  of  thin  films,  f3L  < X,  with  small  tg(s)  there  are  two 
most  important  criteria: 

- influence  of  the  error  of  input  data  on  the  accuracy  of  parameter  reconstruction,  and 

- range  of  parameters,  within  which  accuracy  ofed  by  the  algorithm  is  sufficient. 

To  simulate  experimental  Jnput  data  we  used  (see  also  [3,  4])  a generator  of  random 
numbers  with  normal  distribution: 

Re„  (at,  s)  - 1 (k,  £■)(/  + y j Random) expiargfy1  (k,£)\i  + y 2 Random)), 
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where  - 1 < Random  < 1 . We  studied  the  relative  errors  of  reconstructed  parameters 
and  algorithm  properties  according  to  the  formulas: 


lS^),r  ‘Mm 

L 


%;  A.  = 


'€r  £m  %\  and  A 


I IT 


V -/  ) 

as  functions  of  input  data  errors.  Relevant  curves  are  marked  as  1,2,3  in  Fig. 2, 
respectively.  In  all  experiments  we  accounted  21  frequency  points  in  functional  (2). 

Fig.  2 presents  the  variations  of  relative  errors  when  input  data  errors  for  amplitudes  y, 
and  arguments  y2  change  from  0%  to  10%.  In  Fig.2-a  we  fixed  y2  = 3%  and  in  Fig,2-b 
it  was  y j - 5% . Here  we  had  to  present  curves  2 corresponding  to  the  values  of 
Ac  ■ 100  in  order  to  be  visible  within  the  common  scale.  The  errors  change  randomly, 

however  around  certain  increasing  with  rise  of  errors  level  mean  values.  The  input 
errors  in  arguments  influence  the  accuracy  more  crucially,  and  it  is  clear  that  better  to 
accept  errors  less  than  5%.  Due  to  the  high  accuracy  of  the  algorithm  of  direct  problem 
solution,  there  is  no  restriction  on  reconstructed  parameters  if  one  has  “hypothetical” 
situation  with  accurate  input  data,  i.e.  if  y,  =y:  = 0 . However,  from  numerical 
experiments  we  conclude  that  if  the  error  in  input  data  is  y,  >10%  (for  normalized 
amplitude  that  is  deviation  of  0.1)  and  y2  >5%  (that  is  18°)  it  is  necessary  to  apply 
regularization  to  (2). 
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